Manifestly covariant classical correlation dynamics I. General theory 
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In this series of papers we substantially extend investigations of Israel and Kandrup on nonequi- 
librium statistical mechanics in the framework of special relativity. This is the first one devoted 
to the general mathematical structure. Basing on the action-at-a-distance formalism we obtain a 
single-time Liouville equation. This equation describes the manifestly covariant evolution of the 
distribution function of full classical many-body systems. For such global evolution the Bogoliubov 
functional assumption is justified. In particular, using the Balescu-Wallenborn projection operator 
approach we find that the distribution function of full many-body systems is completely determined 
by the reduced one-body distribution function. A manifestly covariant closed nonlinear equation 
satisfied by the reduced one-body distribution function is rigorously derived. We also discuss ex- 
tensively the generalization to the general relativity especially an application to self-gravitating 
systems. 

PACS numbers: 03.30.-|-p, 52.25.Dg 



I. INTRODUCTION AND OVERVIEW 
A. Some motivations 

The fundamental problem of physics of relativistic clas- 
sical many-body systems is one century old (for a review 
of early investigations see, e.g., Ref. [l|) and, still, re- 
mains one of the most important subjects in studies of 
relativity 0, H, 0, Q ■ At the microscopic level the basis 
of deterministic classical dynamics of relativistic many- 
body systems has been widely explored At the 
macroscopic level relativistic hydrodynamical phenom- 
ena have been well documented [l^, EH- In particular, 
in recent years relativistic effects on macroscopic equi- 
librium and non-equilibrium phenomena have attracted 
considerable theoretical interests, and a variety of re- 
lated subjects are under intensified investigations ranging 
from the internal evolution (or aging) of relativistic sys- 
tems [13, [3 , relativistic generalization of the Maxwell- 
Boltzmann distribution [j, [l3| to relativistic Brownian 
motion [3|. A natural question, thus, is how various 
macroscopic relativistic phenomena emerge from deter- 
ministic relativistic classical many-body dynamics. Long 
time ago it was realized by the Brussels- Austin school 
that, as nonrelativistic many-body systems, the macro- 
scopic evolution generally roots in the creation or an- 
nihilation of correlations between microscopic particles 
as they move in spacetime [H,0- Thus, to investigate 
these issues one has to invent a theory of relativistic clas- 
sical correlation dynamics aiming at understandings of 
the dynamics of distribution function of full many-body 
systems. Unfortunately, this is intellectually challenging 
and remains largely undeveloped. 

From practical viewpoints a theory of relativistic clas- 
sical correlation dynamics may be proven to be power- 
ful in studies of various relativistic transport processes 
in plasma physics JJl and astrophysics tlQj . There 
substantial progresses have been made in past years by 
making use of various relativistic kinetic equations, which 



provide an adequate microscopic description for macro- 
scopic relativistic hydrodynamics. In a pioneering paper 
by Heinz [23] it was realized that the relativistic clas- 
sical kinetic equation may also become a powerful tool 
in studies of quark-gluon plasmas, which has nowadays 
spanned one of the most important branches in high- 
energy physics |2l|. To derive relativistic (classical) ki- 
netic equations there are a number of theoretical propos- 
als. Among of them are two representative approaches: 
The first approach is largely phenomenological, where 
with the help of Stosszahlansatz kinetic equations are for- 
mulated as manifestly covariant versions of their Newto- 
nian analogues so as to suit relativity principles [l^, [ll| . 
The second approach [13, [H, [13, [l^ starts from some rel- 
ativistic BBGKY hierarchy and the cluster expansion is 
further applied. Then, in the crucial step one invokes the 
Bogoliubov functional assumption [1^ to truncate the 
infinite hierarchy. As such, one obtains some relativistic 
kinetic equation. Yet, there does exist the notable excep- 
tion in plasma physics, namely Klimontovich's technique. 
With the help of this technique Klimontovich managed 
to justify the relativistic Landau equation 21] at the full 
microscopic level @, . Since then, it has become 
common to test the validity of various theories by justi- 
fying this kinetic equation [U, [H [S [H [H [H . 

It is fair to say that for some important applications 
such as ultrarelativistic electromagnetic plasmas the par- 
ticle interaction must be described by the quantum field 
theory. Therefore, a complete treatise of relativistic 
statistical effects is required to be built on quantum 
field theories. Then, a general prescription may be to 
generalize-in a manifestly covariant manner-the Green's 
function theory of nonequilibrium quantum statistical 
mechanics, pioneered by Schwinger, Kadanoff and Baym, 
Keldysh and Korcnman [s^, [Ml so as to suit principles of 
quantum field theories. Such a task was undertaken for 
electron-positron plasmas some time ago j32j. The gen- 
eral formalism has been substantially renovated by mak- 
ing use of the functional-integral approach, emerging as 
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relativistic nonequilibrium quantum field theory j33l . |34| . 
These approaches, in combination with the Wigner func- 
tion technique, aUow one to obtain a manifestly covariant 
quantum kinetic equation at the level of the weak cou- 
pling (Born) approximation which, importantly, recovers 
the relativistic Landau equation |27l . l28l | at the classical 
limit. 

Nevertheless there are considerable practical and theo- 
retical reasons which justify to intensify studies of the rel- 
ativistic classical correlation dynamics. The most promi- 
nent one comes the very recent debate on the special 
relativisticgeneralization of the Maxwell-Boltzmann dis- 
tribution triggered by the surprising finding of 
Horwitz and coworkers [1^. In Ref. [2^ it was noticed 
that considerable conceptual and technical differences be- 
tween the relativistic and Newtonian deterministic classi- 
cal dynamics might completely disable simple generaliza- 
tions obtained by the phcnomenological approach. These 
authors proposed a manifestly covariant kinetic equation 
possessing a remarkable mathematical structure differ- 
ent from the one obtained either by other microscopic 
approaches [H, [m, 19, 24,_2&Jpl [H, ^ or by phc- 
nomenological arguments [lOj, . The one century old 
belief-that the Jiittner equilibrium 0, is established 
in relativistic dilute systems-is thereby challenged. At 
this stage it is completely unclear whether such a pecu- 
liar property is specific or universal to relativistic clas- 
sical many-body systems. Thus, it is fundamentally im- 
portant to derive as rigorously as possible the relativistic 
kinetic equation from generic classical many-body dy- 
namics. Then, one may hope to prove or disprove ki- 
netic equations proposed basing on phcnomenological ar- 
guments and, furthermore, to obtain insight into the nu- 
merical experiment on the relativistic equilibrium of clas- 
sical many-body systems Q. Experiences in the New- 
tonian statistical mechanics suggest that the technique 
developed by the Brussels- Austin school [H, [11] may be 
perfectly suitable to fulfill such a task as it manages to 
derive a general kinetic equation [sj without resorting to 
either Stosszahlansatz or the truncation approximation. 
To employ the latter is well-known to be inevitable and, 
in fact, a crucial step in various microscopic approaches. 

There are several important adjacent problems that 
may be explored within the scope of the relativistic clas- 
sical correlation dynamics. One is the relativistic Brown- 
ian motion, which has received much attention in recent 
years (for example, see ^]). Another is to formulate the 
relativistic many-body equilibrium. For this Hakim con- 
jectured long time ago that there exists a Lorentz invari- 
ant equilibrium hierarchy which is invariant under the 
spacetime translation and is merely determined by the 
Jiittner distribution However, to our best knowledge 
there have been no progresses reported so far. 

Finally, new developments of the relativistic classical 
correlation dynamics are urged by the stellar dynamics. 
There, it is generally believed that a satisfactory frame- 
work is provided by the classical (nonequilibrium) statis- 
tical mechanics in the general relativity [18]. In particu- 



lar, the kinetic equation describing the Brownian motion 
of self-gravitating systems suffers from a serious difficulty 
of the infrared divergence [H, [H, [s^ . Such a divergence 
has been explored at the level of the Newtonian physics 
and been known to deeply root in the nonMarkovian and 
collective effect [H, |4l| , the latter of which, interestingly, 
finds the collective dielectric effect in classical plasmas 
[4^ 1 as the analogue. In the Newtonian physics the cor- 
relation dynamics encompasses a standard route towards 
the complete treatise of these effects. It is, therefore, a 
natural hope that the (general) relativistic classical cor- 
relation dynamics may help to understand this issue at 
deeper level. 



B. Nonequilibrium statistical mechanics: 
manifestly and nonmanifestly covariant formalism 

In this series of papers we focus on the (special rel- 
ativistic) classical statistical mechanics. In contrast to 
the relativity in the Newtonian physics the emergence 
of macroscopic irreversibility from deterministic classical 
dynamics nowadays is a well established subject. Indeed, 
for simple chaotic systems especially deterministic dif- 
fusive systems rigorous mathematical investigations, as 
well as various analytical theories [i^ have shown that 
an initial phase space distribution (function) is decom- 
posed into independent components each of which relaxes 
in an irreversible manner (Ruelle-Pollicott resonances), 
and after transient processes the evolution is dominated 
by diffusion at the macroscopic scale. For more compli- 
cated many-body systems (as we switch to in this pa- 
per) a milestone was put by Bogoliubov [2^ who real- 
ized that, there, the evolution is a two-step process: In 
the first step (after transient decays) correlation func- 
tions relax into functionals of one-body distribution func- 
tions. In the second step the one-body distribution func- 
tion evolves following a closed nonlinear kinetic equation, 
and is fully responsible for macroscopic hydrodynamics. 
Later on the concept of the Bogoliubov functional as- 
sumption was justified and substantially extended, espe- 
cially in the profound work by the Brussels- Austin school 
using delicate diagrammatical [s^ and projection opera- 
tor [H, H^l approaches. In particular, the decomposition 
of many-body phase space distribution functions may be 
carried out in the way such that in the absence of in- 
teractions each component corresponds to some degree 
of correlations between particles. These correlations are 
preserved (mixed) in the absence (presence) of interac- 
tions justifying the terminology of "(classical) correlation 
dynamics" . As such, irreversible processes of approach- 
ing equilibrium, as well as the final equilibrium states are 
found to be completely determined by the so-called vac- 
uum of correlations namely the component containing no 
correlations in the absence of interactions. This way the 
kinetic equation is no longer an approximation, rather, 
an exact theory to which the Liouville equation of full 
many-body systems converges in long times. 
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Then, a natural hope is to extend these concepts so 
as to suit relativity principles. At the early stages the 
Brussels- Austin school [l^ li^ undertook systematic at- 
tempts (see Ref. [Ts'l for a review) to generalize the New- 
tonian classical correlation dynamics to special relativity, 
which is built on the canonical formulation of relativis- 
tic dynamics of classical many-body systems Q with the 
Hamiltonian, if necessary, including field degrees of free- 
dom. Unfortunately, as pointed out by many authors 
especially in the notable critical analysis by Hakim 
Israel and Kandrup and Kandrup [l^, to proceed 
along this theoretical line one may have to overcome a 
number of conceptual and technical difficulties. The most 
serious problems are regarding covariance of the theory 
and particularly covariance of the evolution of many-body 
systems. Both deeply root in the absence of absolute 
time in many-body systems. For the former problem it 
should be stressed that by relativity principles it is per- 
fectly legitimate to build any theories in cither manifestly 
or non-manifestly covariant manner. It is exactly the 
latter that the Brussels-Austin school follows. Indeed, 
in the development of old relativistic classical correla- 
tion dynamics a preferred time coordinate is chosen and 
then a Liouville equation is formulated building on the 
Hamiltonian formalism. Furthermore, in an insightful 
work Balescu and Kotera realized that in this framework 
the Lorentz invariance must be understood in terms of 
the Lorentz group action on the distribution functions 
of many-particle phase space (if necessary, enlarged to 
accommodate field degrees of freedom) [13i, i45|]. More 
precisely, the Lorentz group has 10 generators and in the 
group action representation the Liouvillian generates the 
time translation. Then, it is crucial to justify that kinetic 
equations obtained in subsequent calculations (which re- 
quire further delicate approximations) also follow this 
kind of Lorentz invariance. Such a task, if not impos- 
sible, is highly nontrivial and has never been carried out 
so far. For the latter problem relativity principles re- 
quire that the evolution of many-body systems, such as 
approaching equilibrium must be observer- independent. 
Furthermore, it has been a common belief that there ex- 
ists covariant notion of evolution for many-body systems 
P, [H, [H, [2^. Yet, it remains unclear how to reconcile 
these concepts with old relativistic classical correlation 
dynamics. 

A complete different formalism, but equally at the full 
microscopic level, was formulated by Israel and Kandrup 
[Tsl . [igj in the framework of classical general relativity 
(which is of minor importance). The Israel-Kandrup 
formalism, although also starts from deterministic rel- 
ativistic classical dynamics of full interacting many-body 
systems, differs from the old relativistic classical corre- 
lation dynamics in several key aspects. First of all, it 
is manifestly covariant. There, Hakim's analysis of rela- 
tivistic statistical mechanics of /i (namely single particle) 
phase space is substantially extended to F (namely many 
particle) phase space. Then, the underlying relativis- 
tic classical many-body dynamics is formulated in terms 



of the so-called action- at- a- distance formalism (see, e.g., 
Ref. for a review) rather than the Hamiltonian for- 
malism. Remarkably, in such a formalism fields arc not 
considered as independent degrees of freedom. Instead, 
they are carried by complicated retarded (advanced) po- 
tentials. As a result, particles interact in a nonlocal 
manner. And the force acting on given particles must 
be viewed not only as a function of particle's 4-position 
and 4-momentum vector, but also as a functional of the 
world line of all the other particles. Finally, particles are 
treated on different footing in the way that the iV-particle 
system is divided into 1- and {N — l)-particle subsystem 
coupling to each other. As such they managed to em- 
ploy Willis-Picard projection operator approach [4^ to 
obtain a closed nonlinear kinetic equation of one-body 
distribution function. The latter gives the relativistic 
Landau equation [l^, [H, from which, at the full 
microscopic level, Jiittner equilibrium distribution im- 
mediately follows. The success of the Israel-Kandrup 
formalism justified a number of far-reaching concepts. 
Among of them are: the global covariant evolution as 
conjectured by Hakim [l| and the legitimation of building 
relativistic classical non-equilibrium statistical mechan- 
ics on the action-at-a-distance formalism of determinis- 
tic many-body dynamics. The latter nowadays is widely 
accepted [23*, . Importantly, it was suggested that the 
one-body distribution function may, at least asymptoti- 
cally, determine the covariant evolution of the system-the 
very nature of the Bogoliubov functional assumption. 



C. Overview of this series of papers 

Despite of the significant achievements by Israel and 
Kandrup there are many important problems remained 
unsolved. First of all, it is not clear how to go beyond the 
weak coupling approximation within the Israel-Kandrup 
formalism. This issue is of great practical importance 
especially for transport processes in plasmas with elec- 
tromagnetic interactions 01 H^- Also, because of such 
a drawback one fails to predict-at the full microscopic 
level-kinetic equations and (local) equilibrium for dilute 
systems with moderate interaction strength. The solu- 
tion to this problem may solve the puzzle of the rela- 
tivistic generalization of the Maxwell-Boltzmann distri- 
bution which is currently undergoing intensified debate 
B 0, Hi. Then, in Ref. the (local) equihbrium 
was studied under the kinetic approximation. There, 
the principal problem of formulating many-body equi- 
librium 1] has not been attacked. Finally, how to extend 
the Bogoliubov functional assumption rigorously remains 
largely unexplored. 

In this series of papers we substantially extend the in- 
vestigations of Israel and Kandrup [H, [l^ and widely 
explore these issues in the framework of special relativ- 
ity. We present a new, manifestly covariant classical cor- 
relation dynamics. As the basic viewpoint we proceed 
along the line of Refs. [l^ [l^ and formulate the under- 
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lying deterministic classical many-body dynamics in the 
action-at-a-distance formalism with N motion equations 
as fundamental objects [1, 0, Q ■ Such formalism allows 
us to naturally work in physical coordinate and momen- 
tum. The apparent advantage is to admit a theory which 
is manifestly covariant at each step of the manipulations. 
Nevertheless this is by no means merely based on the 
aesthetic viewpoint. Indeed, we have been compelled to 
do so by significant progress achieved recently by both 
physicists and mathematicians, (i) In a 1-dimensional 
numerical simulation it is found that point-like collisions 
tend to drive a relativistic system into equilibrium de- 
scribed by Jiittner distribution Q. At the microscopic 
level the deterministic classical many-body dynamics of 
the underlying system may be described by the action-at- 
a-distance formalism in a rather straightforward manner. 
A statistical theory based on this formalism, in turn, is 
expected to explain the experimental discovery, (ii) The 
Kirpichev-Polyakov theorem partly justifies the long- 
standing conjecture of Hakim on mathematical founda- 
tions of relativistic statistical mechanics It is shown 
that for 1-dimensional relativistic dynamics of classical 
charged systems, which is formulated in terms of the 
Wheeler- Feynmann formalism Q, an ordinary Cauchy 
problem as in the Newtonian mechanics can be stated. 
This is conceptually important because it suggests that 
the global evolution of relativistic many-body systems 
may be formulated in the way analogous to the Newto- 
nian physics, despite that particles interact in a dramat- 
ically different manner. 

The TV motion equations define a natural solu- 
tion space namely the 8A^-dimensional F phase space. 
(Throughout this work the mass-shell constraint is ab- 
sorbed into the distribution functions.) For this phase 
space we may define a distribution function which, re- 
markably, depends on TV proper times. A significant dif- 
ference from the Newtonian physics is that a bundle of 
N world lines, rather than a representation point in the 
F phase space underlies the subsequent analysis of the 
dynamics of distribution functions. Then, we formulate 
TV conservation equations. Moreover, for a large class of 
physical forces such as Lorentz forces the phase space vol- 
ume element is invariant along world lines. As a result, 
despite of the absence of Hamiltonian these conservation 
equations collapse into the manifestly covariant Liouville 
equations. With the TV proper times identified we ob- 
tain a single-time Liouville equation which describes the 
manifestly covariant global evolution of distribution func- 
tion. On this global evolution we may build a manifestly 
covariant theory of classical correlation dynamics by us- 
ing the Balescu-Wallenborn projection operator approach 
[So , [13] . It then follows that the reduced one-body distri- 
bution function fully determines the TV-body distribution 
function, in particular, the entire correlation functions. 
The evolution of the reduced one-body distribution func- 
tion obeys an exact closed kinetic equation. As such, we 
achieve the relativistic Bogoliubov functional assumption 
which, in contrast to earlier theories [H, [UH^, is man- 



ifestly covariant. From the exact closed kinetic equation 
we recover various (Vlasov, Landau, Boltzmann) rela- 
tivistic kinetic equations systematically. The solutions to 
these kinetic equations allow us to pass to macroscopic 
physical observables by carrying out appropriate average 
with respect to them. 

It must be stressed that the present theory cannot 
serve as an alternative to the difficult problem of deter- 
ministic relativistic classical many-body dynamics. Nei- 
ther are they equivalent. As the Newtonian physics in 
order to develop a theory of manifestly covariant classi- 
cal correlation dynamics one actually requires only few 
assumptions regarding dynamical properties of determin- 
istic relativistic many-body systems. These assumptions 
are exactly formulated in this paper. Although to prove 
them may be hard mathematical problems, they have 
apparent physical implications. Thus, these assumptions 
and thereby the present relativistic classical correlation 
dynamics are expected to be applicable for a large class of 
realistic relativistic many-body systems. It is also worth 
pointing out that compared to its Newtonian counter- 
part the classical correlation dynamics presented here 
suffers from additional technical complications. That 
is, the force acting on given particle is determined by 
the world lines of all the other particles and, therefore, 
as one passes from the exact closed kinetic equation to 
special kinetic equations appropriate approximations re- 
garding the world line must be implemented. Indeed, in 
the present work we use the well known relativistic im- 
pulse approximation p^ . [l9j . There, to fully determine 
interactions between two particles at given moment their 
trajectories are considered to be linear. To overcome this 
technical complication one, in principle, needs to either 
expand phase space so as to accommodate particle ac- 
celeration as well as its higher-order derivatives [H, or 
treat fields as independent degrees of freedom which is a 
tractable task. (However, the 1-dimensional relativistic 
dynamics of classical charged systems may be an excep- 
tion, because according to the Kirpichev-Polyakov theo- 
rem the field degrees of freedom are redundant.) Since 
the present work aims at the principle problem of build- 
ing a theory of manifestly covariant classical correlation 
dynamics on the action-at-a-distance formalism, we may 
leave this technical issue at this stage of conceptual de- 
velopment, in particular, if we ignore problems such as 
radiation reaction (for electromagnetic interactions). 

We plan to explore various topics discussed above in 
this and the following [i^ paper. This first one is devoted 
to the general mathematical structure and is organized as 
follows: In Sec.|TT]we first introduce preliminary concepts 
required for developing a theory of manifestly covariant 
global evolution. Then we derive a manifestly covariant 
single-time Liouville equation. In Sec. IIIII the correla- 
tion pattern representation is established in a manifestly 
covariant manner. On this basis we apply the Balescu- 
Wallenborn projection operator approach to the single- 
time Liouville equation. In Sec. IIVI we present the rela- 
tivistic version of the Bogoliubov functional assumption. 
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In particular, we prove the factorization theorem and de- 
rive the exact closed kinetic equation of reduced one-body 
distribution function. We close this paper by outlining 
further applications. In particular, we present an exten- 
sive discussion of generalizing the manifestly covariant 
classical correlation dynamics to the general relativity so 
as to treat self-gravitating systems, which is the main 
subject of Sec. |Vl Some technical details are presented 
in Appendix \X\ and [Bl 

The second paper (denoted as Paper II), written in 
a self-contained manner, is devoted to applications of 
the general principles to relativistic plasmas with elec- 
tromagnetic interactions. The readers interested only in 
applications of the present theory may skip this first one 
and read the second one directly. There, the relativis- 
tic Vlasov, Landau and Boltzmann equation follow sys- 
tematically from the exact manifestly covariant nonlin- 
ear equation which is closed, and satisfied by the reduced 
one-body distribution function. The collision integrals of 
the latter two justify the Jiittner distribution as the spe- 
cial relativistic generalization of the Maxwell-Boltzmann 
distribution. The collective effects, such as the issue of 
correlation at the relativistic many-body equilibrium are 
explored. We summarize this series of work in Paper II. 

Finally we list some of the notations and conventions. 
We choose the unit system with the speed of light c = 1 . 
To distinguish from the Minkowski 4-vector we use the 
bold font to denote the vector in the Euclidean space. 
Greek indices running from to 3 are further used 
to denote the component of the former. The Einstein 
summation convention is applied to these indices. The 
4-dimensional Minkowski space is endowed with the 
metric r]'^" = diag(l, — 1, — 1, — 1) . The scalar product 
of two 4-vectors is defined as a • 6 = r]'^'^af^bi, = a^b'^ . 
In particular, a ■ a = . The arguments of distribution 
functions: {xi,pi,Ti + r) carried by particle i are 
abbreviated as i . In addition to the usual mathematical 
symbols we use the following notations: 



9^, , 



dT,,,. , 



Xi[<;\ , 
di , 



II. 



covariant derivative: = d/dx'^ ; 
volume element in 4-dimensional Minkowski 
space: d^z = dz^ dz^ dz^ dz^ ; 
d-dimensional Dirac function: 
differential form of spacelike 3-surface: 

= ^€^,i,px dx" A dxP A dx^ with e^^pA 
being ±1 when {jjLvpX) is an even (odd) 
permutation of (0123) and being otherwise ; 
world line of particle i ; 
4-position of particle i at proper time ; 
volume element in the /i phase space of 



particle i : di 



Pi 



MANIFESTLY COVARIANT LIOUVILLE 
EQUATIONS 



Let us consider a system consisting of A'' classical point 
particles, each of which has the rest mass m . The parti- 



cles have the action at a distance on other particles. The 
interaction propagates at the speed of light. (To simplify 
discussions throughout this series of papers we shall not 
consider the self-action.) The interactions carry energy 
and momentum to and from particles without the sup- 
port of independent fields. This is the so-called action-at- 
a-distance formalism of relativistic classical many-body 
dynamics pioneered by Fokker ^] and Wheeler and Feyn- 
man [7|. The subject has been widely explored in the 
last half century Q . In this section based on this formal- 
ism we derive a manifestly covariant single-time Liouville 
equation. Analogous to the Newtonian physics it then 
serves as the exact starting point of relativistic classical 
correlation dynamics. 



A. Action-at-a-distance formalism 

The history of the microscopic system is described by a 
bundle of N particle world lines which solve the following 
relativistic motion equations: 



dxl 

dn 

dpi 
dn 



EL 

m 



N 



-F'cxt (a;* , P») + ^ (xj , 



(1) 

(2) 



Here (r^) ,u'^{Ti) ,p^{Ti) are the 4-position, 4-velocity 
and 4-momentum vector of particle i depending on the 
proper time , respectively, F^'^^ is the the external force, 
and F^- is the force acting on particle i by particle j . 
Importantly, we assume that both forces do not depend 
on the acceleration of the acted particle, and consider the 
interacting force F/' with the general form as follows [§] : 



Pl'jiXt,Pz) 



(3) 



-foe 



dTj s{pij ) (o^ij I Pij I ^ij 7 Iji J Cij ) Pvi j 



where T^'^ is an antisymmetric tensor, and the role of 
function s{pij) is to invariantly connect Xi with one (or 
several) points at the world line xj^tj) . The arguments 
in Eq. ([3]) are defined as follows: 



[x, 
(x 



dr-i 



(Xp 



^m) I lij — Pii^M ~ ' 



dp- dp'; 



dTj dTj 



4), 



(4) 



The relation: pi ■ Fij{xi,pi) = is obvious because 
of the skew property of T^^'^ . Such a relation is further 
enforced to external forces: 



Pi ■ Fcy,t{xi,pi) = 0. 



(5) 
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Therefore, the mass-shell constraint: 



(6) 



is guaranteed. To fully determine the interacting force 
Fj'^ relies on the world line of particle j namely Xj(Tj) . 
In Paper II we exemplify the interacting force in the case 
of electromagnetic interactions. 



B. Manifestly covariant single-time Liouville 
equation 

The motion equations (H]) and suggest that we may 
introduce the concept of SA^-dimensional V phase space. 
It is defined as the direct product of N flat p. phase space, 
each of which corresponds to individual particle. That is, 



N 



r 



Mi 



(7) 



where A^^ and Uf are the 4-dimensional Minkowski 
spacetime and momentum space carried by parti- 
cle i , respectively. The volume element of the 
r phase space is dV = OiLi di ,di = d^Xid'^pi . 
There exists a fundamental difference between the F 
phase space defined above and its Newtonian counter- 
part: In contrast to the Newtonian physics Xn = 
(xi(ri),pi(Ti), • • • ,xn{tn),pn{tn)) S T generally does 
not represent a microscopic state because, in general, it 
is not sufficient to uniquely determine the history of the 
full system namely the N particle world lines. Instead, 
in special relativity microscopic states have been conjec- 
tured to be N segments of particle world lines each of 
which is carried by individual particle Furthermore, 
all the microscopic states constitute the space fi^ which 
has the dimensionality > 8N . (Ref. [H] presents a non- 
trivial example where the equality holds.) 

Given a bundle of N particle world lines satisfy- 
ing Eq. ([2]), or more precisely, a microscopic state 
uj € the representation point Xn has the singular 
distribution as YitLi S'^'^H^i - Xi{Ti))6'-'^\pi - Piin)) , 
which depends on ci) through the phase trajectory 
Xi{Ti),pi{Ti) ,i = I,-- - ,iV. By preparing a cloud of 
microscopic states centering at Cj and carrying out the 
average over this ensemble, we expect to smear the 
singular distribution at least along some directions (of 
the F phase space). More precisely, we wish to assert 
the existence of the following 

Assumption 2.1 There exists a probability measure 
on such that the average below 

1^{xi,Pi,Ti, ■ ■ ■ ,xn,pn,tn;xi[(;], ■ ■ ■ ,xn[c;]) 



N 



Y[S^^\x^-x,in))6'-''\p,-p^{n)) 



(8) 



strongly peaks at {xi{ti),pi{ti), ■ ■ ■ , a;Ar(rAr),pAr(rAr)) . 



Remark. Although the details of the average are not 
required for subsequent discussions, it is important to no- 
tice that V possesses N proper times and functionally de- 
pends on the N particle world lines, i.e., , • • • , xat [<j] . 
Physically, it may be interpreted as a probability dis- 
tribution function 49] in the F phase space. Indeed, 
it is obvious that V has the following properties: (i) 
V > and (ii) J dVD — 1 . Given an elementary 
Lebesgue-measurable subset iv = oJi C F , J^dTV 
stands for the probability for particle i to be in LUi in Ti 

Similar to relativistic nonequilibrium mechanics built 
on the /J, phase space [ij, from the causality follows 



Lemma 2.2. The distribution function T> satisfies 
lim I? = 0, Vi 

Ti — *-±00 

in the sense of Lebesgue measure. 



(9) 



Proof. Consider the probability density at the proper 
times (ti, • • • , tat) defined as 



N 



p{xi,Ti, ■ ■ ■ ,xn,tn;xi[(;], ■ ■ ■ ,xn[';]) = jY^d^PiVl^lO) 



Then, the probability for particle i to stay inside a finite 
Lebesgue-measurable volume oJi € Aij at the proper time 
n {i = 1, - ■ ■ ,N) is 



d'^XNP- (11) 



n{Ti, ■■■ ,tn) = d xi 



Notice that for given i the world line Xi [<;] , having crossed 
uji , could not return. Thus, 



lini n(ri, • • • ,tn) ^0 

— >±00 



(12) 



for arbitrary LUi . From this the lemma follows. Q.E.D. 

However, the distribution function T> has a remarkable 
property. That is, with the momenta integrated out it is 
reduced into a distribution function defined on the entire 
Minkowski spacetime. Thus, the distribution function T> 
is not physical. Instead, to match macroscopic obser- 
vations we have to further introduce the physical distri- 
bution function Af such that with momenta integrated 
out it is reduced into a distribution function defined on 
a spacelike 3-surface-a subspace of the Minkowski space- 
time. Mathematically, it is defined as 

N 

N{xi,pi,-- - ,XN:PN\xi[(i], - ■ ■ ,xn[(,])= j'^dnV ili) 

•' i=l 

which is the straightforward generalization of its p, phase 
space counterpart [l|. The normalization of J\f follows 
from 
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Lemma 2.3. The physical distribution function J\f 
satisfies 



/ dT,^^d'^piu1; ■ ■ ■ dY.^.^d'^pNU^Af 



1 , 



(14) 



Here Si is an arbitrary spacelike i- surface in Aif with 
the vectorial value of the differential form as dS^. . 

Proof. For the volume element d'^Xi we have 

dx'^ 

d*Xi ^ dT,fj,^dxf ^ dTidT,f^^—-^ , (15) 

where dS^^i is the differential form of an arbitrary space- 
like 3-surface S,; : 



dE,, 



1 
3! 



e^i^pA dXf A dx'^ A dx^ 



(16) 



Here e^^p\ is ±1 when {^vpX) is an even (odd) permu- 
tation of (0123) and is otherwise. Then, 

d^x.d^p^S^^^x, - x,{n))S^^'>{p,-p,{n)) (17) 
= dndJ:^^^d^p,S^^Hx^ - xdri))6^^\p, - p^in)) 

an 

= dTidY.fj_.d'^Pi u'^ d''^\xi - Xi{n))S'^'^'> {pi - Piin)) . 



Notice that J dTV = 1. We then substitute Eq. Hit 
into it and integrate out r^'s. With Eqs. ([5]) and (fT51 
taken into account eventually we obtain Eq. (I14p . Q.E.D 



Remark. This lemma shows that the physical distri- 
bution function Af is the analogue of the F phase space 
distribution function in the Newtonian physics. The 7N- 
dimensional manifold on which Af is normalized may be 
considered as the effective phase space. For the conve- 
nience below here we also introduce the space of physical 
one-body distribution Ti. defined as: 



n 



h\h : ^ (g) U"^ ^ R+ U {0} , 



lim N^^ 

Af— » + oo 



dY.^d^pu^'h^l} . (18) 



To further proceed we introduce the Liouvillian £ 
which is decomposed into the free and interacting part, 
i.e., £ = £P+X£,' with £P and A£' the free and interacting 
Liouvillian, respectively. They are defined as 



£" = -E 

1=1 



d_ 



,(19) 



where XC[j is the two-body interacting Liouvillian. 
Notice that here the dimensionless parameter A charac- 
terizes the interaction strength. Now we are ready to 
prove the following Liouville theorem which justifies the 
covariant notion of evolution in the F phase space, and 
constitutes the exact starting point of the succeeding 
sections: 

Theorem 2.4. If both the external and interaction force 
are conservative, i.e., 

—F:^,{x,,p,) = 0, —Ftjix^p,) =0, (20) 

then the distribution function 'D{xi,pi,Ti + 
T,- ■ ■ ,xn,pn,tn + t;xi[(;],- ■ ■ ,xm[<^]) satisfies the 
following Liouville equation: 



d 



-£ P = 0. 



(21) 



Proof. Suppose that at the proper times 
(ti,''' i^n) particles are in the volume lu — 
d'^Xi{(;i)d'^pi{(;i) with a probability density 
I?(a;i,pi,Ti -I- T, • • • ,XN,PN,TN + t; Xi[(i], ■ ■ ■ ,XAr[^]). 
Now let some particle say i evolves from q to along 
the world line Xi[(;;] , while others frozen. As such the 
probability is conserved. That is, 

d'^Xi{(;,)d'^Pi{(;,)Y[d'^Xjd'^Pj'D{- ■ ■ Xi{<;i),pi{(ii),Q ■■■) 

= d'^Xi{<,'i)d'^Pi{(;'i) Jl d'^Xjd'^pjV{- ■ ■ Xi{c;l),p,{4),c;l ■■■) 

(22) 

with all the irrelevant arguments in V suppressed. Since 
the forces are conservative, the Lebesgue measure is con- 
served, i.e., d'^Xi{<;i)d'^Pt{<;i) = d'^Xt{^-)d'^pi{<;-) . Eq. 
then gives 



■ ■ a;i(Q),K(<r»),Q • • • ) = X'(- • • Xi{4),Pt{4),4 • ■ ■ ) 



(23) 



Letting — s- we obtain 



N 



{xt.Pt) + ^Fl'.{x,,pi) 
I?(a;i,pi,^i, • • • ,XN,PN,'iN;xi[i,],- • • ,a;Ar[<r]) = . 



d 



(24) 

Totally we have N such Liouville equations with different 
evolution time q . 

Furthermore, we let the N proper times change uni- 
formly, i.e.. 



dc;i 



di,N — dr 



(25) 
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or equivalently, q = + r with t^'s the initial proper 
times. Then, adding the TV LiouviUe equations together 
Eq. jni) foUows. Q.E.D. 

Remarks, (i) In the absence of interactions the above 
single-time LiouviUe equation was derived by Hakim 
The manifestly covariant LiouviUe equation originates at 
the probability conservation law, and is irrespective of 
the Hamiltonian character of relativistic classical many- 
body dynamics. Of course, in the Hamiltonian formal- 
ism the (single-time) LiouviUe equation trivially exists 
albeit nonmanifestly covariant [l5l l45j . In the action-at- 
a-distance formalism the classical many-body dynamics 
has non-Hamiltonian character. To justify the LiouviUe 
equation the constraint of conservative force must be im- 
posed. Provided that such constraint is released one may 
enlarge the F phase space and subsequently arrive at a 
generalized single-time LiouviUe equation 

(ii) The single-time LiouviUe equation shows that the 
distribution function V depends on a proper time t 
parametrizing the global evolution in the F phase space. 
Yet, it still depends on the other — 1 proper times 
associated with the initial condition of the global evo- 
lution. Upon passing from N LiouviUe equations ((24l) 
to the single-time LiouviUe equation (pij) some details 
of deterministic classical many-body dynamics are lost. 
Thus, the single-time LiouviUe equation is not equivalent 
to relativistic motion equations. The peculiar feature of 
multiple proper times roots in the noninstantaneous fea- 
ture of relativistic force and causes a conceptual differ- 
ence from the Newtonian physics. In the latter case the 
instantaneous nature of forces allows one to parametrize 
the global evolution by the usual coordinate time. Conse- 
quently, in the sense of dynamics of distribution functions 
(single-time) LiouviUe equation is equivalent to motion 
equations. 

(iii) For Eq. (|2ip let us integrate out r^'s. Taking into 
account Eq. we find 

£7V(xi,pi, • • • ,xn,pn;xi[^], ■ ■ ■ ,xn[<,]) = 0. (26) 

Provided that the functional dependence of Af on 
xi[<^], • • • ,a;Ar[<;] is released it then recovers the exact 
starting point of Refs. llM [l9l (in the framework of gen- 
eral relativity) and Ref. [2J]. The very nature of such 
functional independence on the N particle world lines 
is that at given global proper time t the representa- 
tion point is sufficient to uniquely determine the particle 
world lines. This does happen in the asymptotic sense. 
It turns out that, as the Newtonian physics, for suffi- 
ciently large proper time r the relativistic many-body 
systems may also lose the memory on "initial" world line 
segments-the profound change of deterministic classical 
dynamics. Consequently, the instantaneous evolution of 
A/" is determined by specified phase trajectories associ- 
ated with N particles each of which, importantly, merely 
depends on the phase coordinate of given particle. In 
Sec. II VI we will show that the single-time LiouviUe equa- 
tion (I21D. indeed, admits such kind of solution in the 



thermodynamic limit: N —f +oo . However, we are not 
able to justify this picture for finite but sufficiently large 
N. 



III. CORRELATION DYNAMICS ANALYSIS OF 
GLOBAL EVOLUTION 

In this section we wish to apply the Balescu- 
Wallenborn projection operator approach [H, H^l to the 
single-time LiouviUe equation (j2ip . In doing so we 
achieve a theory of relativistic classical correlation dy- 
namics in the way that at each step the manipulations 
are manifestly covariant. With this approach we man- 
age to split the proper-time dependent distribution func- 
tion into the kinetic and non-kinetic component. Each 
of them independently evolves, and the latter decays for 
sufficiently large global proper time. 

It should be emphasized that given appropriate as- 
sumptions regarding the deterministic classical many- 
body dynamics, as to be specified below, (In other words, 
we will deal with some axiom dynamical system.) all the 
results presented in Sec. IIIII and IIVI are mathematically 
rigorous. 



A. Reduced distribution function representation 

Let us first introduce the concept of reduced distribu- 
tion function. From now on we denote the arguments: 
{xi,pi, Ti + t) of distribution functions as i . Then, the 
reduced s-body distribution function is defined as 

Vsiii, ■ ■ ■ ,is;a;.,J<r], • • • ,a;ij^]) 

N 



J = S + 1 



yi<ii<---<is<N. (27) 



Notice that here the particle groups (ii,--- ,is) and 
(is+i,--- jiw) constitute a partition of the full system 
(1, • • • , N) . Moreover, the reduced s-body distribution 
function Vg stands for the probability density for particle 
ij to be at (xi^ , pi^ ) at the proper time Ti. +t (1 < j < s) . 
The normalization can be shown to be 

/ Y^dij'Dsiiir ■ ■ , i^; a^ij^], • • • ,Xi^[<;]) = 1. (28) 

Then, for Eq. ((2T|) we integrate out the phase coordi- 
nates of particle is+i, ■ ■ ■ , «Ar • As a result we obtain 



d 



J = l 3<j' ) 

^ J dis+i \C',^,^^Vg+i ,(29) 



is+i — 1 



is+i / ii 
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which is a manifestly covariant relativistic BBGKY hi- 
erarchy. Here in order to make the formula com- 
pact we have omitted all the arguments of Vg and 
■ Notice that Eq. differs from that studied 
in Refs. [13, [H,!!^, where nonmanifestly covariant rel- 
ativistic BBGKY hierarchies were derived for physical 
distribution functions. 

Let us further define the distribution vector: 

3 =({I?i},p2},--- ,PA'} = 2?). (30) 

Notice that, for each {Vg} there are N\/[{N — sljsl] com- 
ponents each of which corresponds to a s-particle group 
(ii, • • • , is) with 1 < ii < ■ ■ ■ < is < N . With the help 
of this definition the BBGICY hierarchy is rewritten in a 
compact form: 

(|:-£)i(.o, (31) 

The projection to general s-particle states, denoted as 
(«i , • • • , is I is identical to Eq. provided that the ma- 
trix elements are set to be 



k=l 




and to be zero otherwise. Notice that the two-body Liou- 
villian XC, , is determined by the world lines of parti- 
cles ifc and Zj+i . Eqs. may be considered as the 

reduced distribution function representation of the single- 
time Liouville equation pT|) . 

Likewise, we may carry out the same program for the 
physical distribution function A/" and cast Liouville equa- 
tion ([^5]) into another manifestly covariant BBGKY hi- 
erarchy. The details are presented in Appendix VK\ 



B. Correlation pattern representation 

To proceed further we introduce the so-called cor- 
relation pattern representation. Consider a general s- 
particle object F{1, ■ ■ ■ , s; Xk^ M, • • • , [i]) possessing 
the permutation symmetry with respect to 12 • • • s , where 
1 < ki < N , and ki ki' ,y i ^ i' . A unique correlation 
pattern, denoted as jr^) (or {Ts\) may be assigned to it 
in the following way: If F is factorized into s components 
each of which merely depends on i = {xi,pi,Ti + r) and 
XkA<^] {1 < I < s), then |r,) = \0s) = |l|2|---|s) (or 
{Ts\ — (l|2|---|s|), which is called vacuum state . In 
general, F is at most partially factorized and called cor- 
relation state. And the correlation pattern describes the 
factorization structure. More precisely, suppose that F 
possesses the factorization as follows: 



j 

= n ^^i^i'^k,. H , Xk,^+i H , • • ■ , a;fe,.+,. H) , (33) 

i=l 

where Pi = (s^, -I- 1, • • • , -f ji) 7^ , i = 1, 2, • • • J < 
N constitute a partition of (1, 2, • • • , s) : 

PiU---UPj = (I,-- - ,s), P,nP,. =0,Vi ^i'. (34) 

Then, the correlation pattern is jPs) = |Pi IP2I • • • |Pj) (or 
(r^l = (P1IP2I • • • |Pj|)- Notice that the correlation pat- 
tern possesses the permutation symmetry with respect 
to P1P2 • • • Pj • Therefore, two correlation patterns are 
considered to be identical if they differ only in the order 
of P,. 

With this definition the distribution vector may be de- 
composed in a more delicate manner: 



Vsiii, ■ ■ ■ ,Zs;XiJ^], • • • ,XiA'i]) 
= X! \^s){Ts\T>s{ii, - ■ ■ ,is;xt,['^],- ■ ■ ,XiJ^]) 

s 

= E E |Pi|P2|---|Pj)(Pi|P2|---|Pj|i?.(zi,--- ,^s;xnH,--- ,^..H) (35) 
j=i Pi-Pj 



for general reduced s-body distribution function, which 
is the reformulation of the cluster expansion. Here 
Pi , • • • , Pj is a partition of (ii • • • i^) . Notice that in the 
third line any two terms do not possess identical corre- 



lation pattern. Furthermore, in Appendix [B] we give the 
matrix elements of the preliminary operators and A£' 
in this representation. It should be stressed that such 
a cluster expansion differs from traditional one (36j in 
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that the distribution functions depend on particle world 
lines. Such a concept was first introduced by Klimon- 
tovich in the Newtonian context [l^ and was generalized 
to special relativity-in a manifestly covariant manner-by 
Hakim [li]. 

In this way we have assigned all the preliminary quan- 
tities to be used below, i.e., S) , £P and A£' a unique 
decomposition in the correlation pattern representation. 
Thus, we may define the following vacuum and correla- 
tion operator, denoted as V and C , respectively: 



vir. 



So. 



(l-Jo.rJ|r,) (36) 



which are diagonal in the correlation pattern representa- 
tion jr^) (and also in (Frl). It is easy to show that V 
and C constitute an orthogonal decomposition, i.e., 

V + C = l, 
= V,C2 = C,VC = CV = 0. (37) 

With the correlation pattern representations of 
D , £^ A£' , V and C as building blocks, one may 
proceed to establish representations of more complicated 
operators with the help of following properties: 

1. Completeness. J2r^r |rr)(rr| = l; 

2. Orthogonality. {Tr\T'^,) ^ Srr' Sr.r' ! 

3. Linearity. For arbitrary operators A and B , {A + 

B)\rr) ^ A\rr) + BlTr) . 



C. Propagating operator and irreducible evolution 
operator 



From now on we consider closed systems namely F^f^^ — 



. In the absence of interactions, Eq. (|2T|) is reduced into 



I- - £" ] I? = . 

OT 



(38) 



In this case, the evolution is determined by the propagat- 
ing operator il°(T) , or equivalently, its resolvent Dl'^{z) 
which are defined as 

£t"(r)=e^^\ ^0(z) = — (39) 

They are related through 

lt°(T) = (2^)-i / dze-^""fR"(z), (40) 
Jc 

where the contour C lies above all the singularities of the 
Laplace transform of V . 

We further introduce the irreducible evolution operator 
€(t) and its Laplace transform £{z) : 

oo 

£{z) = ^A"+i£'{C?l"(z)£'}", (41) 

£(r) = (27r)-i / dze-'^^£iz). 
Jc 



Then, the axiom dynamical system that we will consider 
throughout this work is defined as such systems that 
satisfy the following 



Assumption 3.1. The operators: Vf(z)V, 

Y£{z)m"{z)C, Cm°{z)£{z)Y and Cm°{z) + 
C$K'^(z)£(z)9^*'(z)C are regular and nonvanishing 
at z — . 

We stress that this assumption merely introduces restric- 
tions on microscopic interactions. A heuristic example 
where this assumption is applicable is a system com- 
posed of identical particles which interact with each other 
through short-ranged interactions. 



D. Evolution of the kinetic component of 
distribution vector 



Under Assumption 3.1. the following two theorems 
are straightforward generalizations of their Newtonian 
counterparts, and we shall not present the proof here [11]. 

— ^ 

Theorem 3.2. The distribution vector S may be de- 
composed into the kinetic component Ilk. 2) and the nonk- 
inetic component lin k. 2) , *-e-, 



(42) 



The vacuum part of the former satisfies a closed hierar- 
chy: 



d_ 



- VrV ) VHk.S =0, 



(43) 



where the operator VFV is given by the following func- 
tional equation: 

VFV = V£V+ / dsVg(s)Vexp(-sVFV), 
Jo 

Vq{s)V = (27r)-i / dze-'"'Y£iz)m"iz)C\£'V . 
Jc 



(44) 



Remark. The operators Hk. and Hn.k. are r-independent 
nonlinear functionals of the interacting Liouvillian A£' . 
However, their explicit expressions are not needed for 
subsequent analysis and, therefore, we shall not present 
them here. Then, carrying out the Laplace transform 
[denoting the Laplace transform of g(s) as §2] and sub- 
stituting Eq. dH]) into Eq. we find 



Vfl^V = ^ A"+iV£'[CK"(z)£']"V. (45) 

n=l 
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Notice that the operator VTV depends on the given par- 
ticle world lines through the interacting Liouvillian A£' . 

For the correlation part, as shown in the following 
theorem, it is fully determined by the vacuum part: 

Theorem 3.3. The kinetic component of the distribution 
vector Ilk. S satisfies 

Cflk.S) = / ds rfs'Clt°(s-s')e(s') 
Jo Jo 

X exp(-sVrV)nk. S . (46) 



Eqs. and constitute the main equations of the 
manifestly covariant classical correlation dynamics. Cru- 
cially, they differ from traditional correlation dynamics 
[Til . UM . [361 1 in that both the operators T , € , Ilk. and the 
distribution vector S) are functionals of particle world 
lines. Such significant difference causes additional con- 
ceptual and technical complications. 

IV. BOGOLIUBOV FUNCTIONAL 
ASSUMPTION 

In the thermodynamic limit N +oo the kinetic com- 
ponent obeys an infinite hierarchy. In this section we 
show that such an infinite hierarchy may be reduced into 
a closed highly nonlinear equation which, remarkably, ad- 
mits a solution uniquely determining the physical distri- 
bution function of full many-body systems. In doing so 
we achieve a manifestly covariant version of the Bogoli- 
ubov functional assumption. 

A. Hakim-Israel-Kandrup trial solution 

For the convenience below we introduce the mapping 
X which uniquely converts a representation point in fi 
space into a world line. It is defined as follows: 

X : e A* i-> -^(x.p) = x[s] (47) 

with the world line x[s] satisfying 

(i) a;|,=o ^ x, 

(ii) v., v%5<l; 

(iii) d^xd^p is invariant along the phase trajectory: 
{x[s],m£x[s]) . 

The property (i) implies that x[s] is a 1-dimensional 
submanifold of which passes through x with the 
4-momentum vector p , (ii) implies that the world line 
^(x.p) preserves the causality, and (iii) implies that par- 
ticles move in the way as exposed to fictitious external 
field which introduces "conservative forces" . Notice that 



the details of X , although complicated generally, are 
unimportant at this stage. 

Eq. (|43| is a closed hierarchy of Ilk. S . Thus, the latter 
as a whole may be regarded as a new distribution vec- 
tor, denoted as J)" = Ilk. S , which possesses the similar 
structure as S , i.e., S° = ({P?} , ,•••). In this 
part we come to study a particular solution of 2)'^ with 
the general component as 

V°{ii, ■ ■ ■ ,is;x,,[<;], ■ ■ ■ ,a:,J<?]) 
= 2?°(ii,--- ,is; A'(2;.^_p.j, • • • ,A'(^.^^p.j) 
= pO(ii,... ,A',J. (48) 

Also, we define the distribution vector J)oo as follows: 

7v^"(r,(zi,.- - ,z,)|s°^(r,(zi,... ,z,)|Soo, 

v(r,(zi,... ,z,)|,J>l. (49) 

In Ref. [l[-in an implicit manner-Hakim noticed that 
in order to derive kinetic equations the distribution func- 
tions in the 8-dimensional /i phase space have to func- 
tionally depend on the particle world line, and the lat- 
ter merely relies on the phase coordinates in the sin- 
gle particle m phase space. In the work of Israel and 
Kandrup [la, this concept was reinforced and ex- 
plicitly formulated. Such an idea lies at the heart of 
Eq. (|48|) . For this reason we may call TIq with the 
components given by Eq. (|48p the Hakim-Israel-Kandrup 
(HIK) trial solution. From the HIK trial solution an im- 
portant fact follows: In general, to determine expres- 
sion of Eq. (|43|l explicitly requires calculations of the 
general matrix element • • ■ exp{TVrV}|ii | • • • 1?/) 
with j < j' , where the sequence ii-'-ijorj' satisfies 
is is' ,y s ^ s' . On one hand, exp{rVrV} explicitly 
depends on the particle world lines given by the acted 
distribution vector. On the other hand, the particles, 
though having identical phase coordinates, may still be 
"distinguished" through their world lines which generally 
are coupled to each other. Consequently, the quantity: 
{ii\---\ij\exp{TYTV}\ii\---\ij,){ii\---\iy\^" is sensi- 
tive to particles Zj+i ■ ■ ■ ij in the intermediate states. The 
substantial simplification introduced by the HIK trial so- 
lution is just to remove this sensitivity by decoupling the 
particle world lines. That is, 

(zi|---|i,|exp{TVrV}|*i|...K,|*,+i|...|z,.) 

= (*i|-..|*,|exp{Tvrv}|*i|...K,|*;.+i|...|z;.,) 

x(*i|---n,n;+iI---ivis°, (50) 

where the particle groups (ij+i,--- ,ij>) and 
{i'j+i, ■ ■ ■ li'ji) are not identical. 

For the HIK trial solution the factorization theorem 
below shows that Eq. (j43]) may be reduced into a 
single closed equation of reduced one-body distribution 
function, which is the generalization of Clavin's theorem 
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in the Newtonian physics [50|. Furthermore, the stationary solution with respect to the 

T -parametrized evolution above may be reduced into the 

Proposition 3.4. In the limit N — s- +00 the infinite following closed equation: 
hierarchy J^gp may be reduced into a single closed 
equation of reduced one-body distribution function as 

^"^^""""^ - / d^,,d%ul^XC[,f{2)] /(I) = K[f] 

Q^+<d^i- J d2\C[,Vi2;X2)jVil;X,) (52) 

y"'^2---y"dj(l|V(r-£)V|l|---|j) J]^2?(s;A',). ^^^^ f e H a physical one-body distribution and K a 

s=i nonlinear functional of f . 

(51) 

I 



Proof. Let us consider the HIK trial solution. We may formally solve the hierarchy (|43|) and project the solution 
D'^(t) to the one-body vacuum state say {i\ . Consequently, we obtain (In order to make the formula compact for the 
moment we suppress all the arguments of except the proper time r parametrizing the global evolution.) 

N N 

(z|SO(r) = (*| exp{rV^V}|^)(^|^°(0) + E E exp{rV^V}|^|^l | • • • \^,){^\^,\ ■ ■ ■ |»,|S"(0) , (53) 

j=l 

where the sequence iii - ■ ■ ij satisfies is ^ i ,y s and is ^ is' ,Vs s' . Without loss of generality the sequence is 
ordered in the way that it starts from i . Because of the particle symmetry namely Eq. (|50p we may simplify Eq. (j53p 
as (Without loss of generality we set i = 1 .) 

(l|SO(r) = E I^TTyl exp{rVrV}|l|2| • • • b-)(l|2| • • • ^"{0) , (54) 

where the combinatorial factors arises from the particle symmetry. Notice that the matrix element of the operator 
expjrVrV} is well defined in the limit N +00 because the world lines are given by the mapping X . Multiplying 
both sides by N and taking Eq. (|49l) into account we find 



(l|S„o(r) = E exp{rVrV}|l|2| • • • |j) (1|2| • • • |j|S.o(0) (55) 

i=i 

from Eq. 

In general, for finite N projecting S)°(t) to a vacuum state, say (zi| • • • , 1 < j < N we obtain 
(zi|---|z,|S"(t) = (^l|•••|^,|exp{rV^V}V^"(0) 

N N 

= E E (H|---N,|exp{rVrV}|ji|...|z,|z,+i|...|*„)(zi|...|z,|z,+i|...|z„|SO(0), (56) 

where for fixed n the sequence ii---in satisfies is ^ is',Vs =/= s' , 1 < s,s' < n. In the second equality for the 
intermediate vacuum state the particle order again needs to be distinguished and, without loss of generality, we set 
the leading j particles to be iii2 ■ ■ ■ ij ■ Because of the particle symmetry with the limit N +00 taken we obtain 
(setting ifc = fc,A;= I,-- - ,j) 

(iI---IjISoo(t) 

00 

= E(l|---b1exp{rVrV}|l|...|n)(l|...|n|S^(0)^i? (57) 

n=j 

or equivalently, 

Soo(t) =exp{TVrV}So,(0) (58) 
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by noticing (1| • • • |j| exp{rVrV}|l| • • • In) =0 for j > n . 

We then come to analyze R . For this purpose we consider the diagrammatical presentation of general matrix 
element: (1| • • • \j\ exp{TVrV}|l| • • • \n) , j < n describing the transition from the intermediate to final vacuum state. 
Fully parallel to the Newtonian physics 36] , starting from Eq. (|^^ and employing the correlation pattern representa- 
tion of operators il*^ and X£,' (see Appendix [B| one may show that diagrams representing the transition element are 
composed of j disconnected clusters. Each of them involves a particle group Pj (i = 1, 2, • • • , j) with 1 + ji particles 
(ji ^ 0) which are labeled as i and Si, Si + 1, ■ ■ ■ , Si + ji , j + 1 = si < ■ ■ ■ < Sj < n . Pj's constitute a partition of 
(I,... ,n): 



PlU---UPj = (I,--- ,7l), ieP^ 

p^nP,' = 0, yi^i'. 



(59) 



These clusters share a common diagrammatical feature: They start from 1 +ji disconnected particle lines at the right- 
most side, which gives the vacuum state |Oi(Pi)) = |j|si|si + l \ ■ ■ ■ \si + ji) . Upon propagating to the left particles 
Si, Si + 1, ■ ■ ■ , Si+ ji are annihilated and eventually at the left-most side only particle i is left giving the vacuum state 
{i\ . Thus, we factorize the matrix element R into j components: 



n=i Pi---Pj 1=1 



= E ■ • • E n exp{rVrV}Ki| . . . |z,.)(*i| • • • |z.JS.o(0) 

Si — 1 Sj —1 i—1 

j oo 

= HE (*|e^p{^vrv}|*i|...|z,j(*i|-..|z,jSo,(o), 



(60) 



= 1 Si=l 



where in the second equality we use the particle symmetry namely Eq. (|50p to make the change of variables: P^ 
*2, ■ ' ■ J *si in such a way that ii = i and ig ^ is' , s ^ s' . 
With Eq. dini) substituted Eq. dST]), together with Eq. dSS]) then gives 



3 oo 



(1| • • • |j|2)oo(t) = n E exp{rVrV}|zi| • • • 

i=l Si = l 

I 



|Soo(0): 



V j > 1 



(61) 



Such an infinite hierarchy is solved by 

(l|...|j|Soo(r) = n Ari??(s;A',) 

'j 

= Y{'b{Xs,Ps,Ts+T-X(,^,p^^), Vi> 1,(62) 



where in the last equality we have retrieved all the ar- 
guments of the distribution functions. Indeed, with such 
a solution inserted the infinite hierarchy (|6ip is reduced 
into a single closed equation as follows: 



V[l;Xi) (63) 

oo j 

J2 (1| exp{TVrV}|l| • • • \j) f] V{s; Xs)\r=o ■ 

j=l s=l 



Taking the derivative with respect to r wc obtain 

— P(l; A-i) = V (l|VrVexp{TVrV}3oo(0) 



^ (l|VrV|l| • • • b-)(l| • • • IjI exp{rVrV}S„o(0) 



5](i|vrv|i|...b-)n^(^;'^^)' 



(64) 

which is the differential form of Eq. ([55]) . We thus prove 
the first part of the proposition. 

Let us now fix a spacelike 3-surface Si for particle i 
and denote the coordinate as Oi . Then, the stationary 
I?(i; A'i)-with respect to the r-parametrized evolution- 
assumes the following form: 



Vi^Xi) = f{ai,pi;X(^^^^p^^{Ti - <;i)) . 



(65) 



Here q is the proper time when the world line 
'^ixi,pi) passes through the phase point {xi,pi), and 



14 



/{(^ijPi'i ^{xi,pi)is)) is a distribution function peaking 
at the spacelike 3-surface . Inserting Eq. (p5)) into 
Eq. ([5T|) and integrating out q , i > 2 and n , we obtain 
a general kinetic equation ([5^ with the collision integral 
given by 



K[f] = / dT,f^2d'^P2uf^ ■ ■ ■ dT,^,jd*pju'^ 



j>2 -i^^mt 

x(l|V(r-£)V|l|...b-) n/(*), (66) 

i=l 

where 

f{xi,pi)=J ds f{ai,pi;X(a:i,pi){s)) eHi. (67) 

Notice that in Eq. ((66| the two-body interacting Liouvil- 
lian XC[j{xi,pi;Xj,pj) is a functional of the world lines 
p.) and X^^. p.^ . The second part of the proposition 
then follows. Q.E.D. 

Remark, (i) By definition the distribution function 
([65]) satisfies 2?|ri^±oo = and thereby is compati- 
ble with Lemma 2.2. / may take a particular form: 
fxixiyPi) 6{Ti - <;i) with fKixi,Pi) e "Hi, which was 
given-in an implicit manner-in Ref. [l9!|. (ii) An equa- 
tion similar to Eq. (jSTjl has been obtained by Hakim [l| 
for dilute gases with scalar and electromagnetic interac- 
tions. There, the equation is derived under the weak 
coupling approximation (namely the second order inter- 
action expansion). In contrast, provided that the force 
is conservative Eq. (jSip presented here is exact including 
all the higher order (three-, four-body, etc.) correlations. 
It encompasses a route to go beyond the weak coupling 
approximation for dilute gases, and is essential to the 
justification of the relativistic Boltzmann equation (to 
be detailed in Paper II). 



B. Physical correlation functions 

The solution of Eq. (|52p fully determines the entire 
physical correlation functions. Indeed, let us replace 
Ilk. 2) in Eq. with D ° , the components of which 
are given by Eq. (|48| . Taking into account Eq. ((65)) . by 
straightforward calculations wc may find the hierarchy of 
physical correlations in the limit: N +oo , denoted as 
COloo , to be 

Cmoo = / ds ds'CU"{s- s')^{s')V 
Jo Jo 

X exp(-sVrV)^oc. ■ (68) 

In the derivation we notice that the vacuum part of the 
physical distribution vector is given by 



(69) 



for arbitrary vacuum state • • • |jn| • 

Alternatively, we may project Eq. (j68p to the general j- 
particle correlation pattern Tj{ii, ■ ■ ■ ^ • • • , ij) 

and arrive at 

{tSoo^Y. / '^^ ds'llfit,) (70) 
x{{rj\CLi°{s - s')e(s')Vexp(-sVrV)|ii| • • • , 

where (•) stands for the replacement: J di — s- 
J^^^jji d^^id'^Piu'^ for the annihilation vertex [see 

Eq. (IMI) ]. 

This exact formula provides a principle to calculate ar- 
bitrary physical correlation functions, provided that the 
physical one-body distribution function namely f{x,p) 
is given. Eqs. ((52)) and ([70]) show that (after transient 
processes) physical correlation functions relax into func- 
tionals of physical one-body distribution function. The 
latter obeys closed kinetic equation ([5^ from which a rel- 
ativistic hydrodynamic description stems. Thus, we jus- 
tify the manifestly covariant the Bogoliubov functional 
assumption. 



SOME REMARKS ON FURTHER 
APPLICATIONS 



It has been a long standing problem to reconcile sta- 
tistical mechanics and relativity principles. In recent 
years this subject has become fundamentally important 
to studies in many fields as mentioned in the introductory 
section. In 1984 Israel and Kandrup made substantial 
progresses in this direction [l^, [l^. There, the authors 
formulated manifestly covariant classical nonequilibrium 
statistical mechanics and successfully applied it to rel- 
ativistic plasmas with electromagnetic interactions. In 
particular, the relativistic Landau equation 27[ was jus- 
tified at the full microscopic level. It differed remark- 
ably from earlier attempts [l^ [2^ in that the theory is 
manifestly covariant at each step of the manipulation. 
As such, the formulated nonequilibrium statistical me- 
chanics (particularly various approximations unavoided 
on top of it) is guaranteed to suit relativity principles au- 
tomatically. That impeding further development lies, as 
pointed out in the original paper (l9j , in that the Israel- 
Kandrup formalism fails to go beyond the weak coupling 
approximation and to capture the collective effects. (It is 
important that the latter heals the well-known logarith- 
mic divergence of the Landau equation.) One of the main 
purposes of this series of papers, indeed, is to attack this 
difhcult problem. 

In this paper we start from the single-time Liouville 
equation (PT|) and derive rigorously the manifestly co- 
variant closed nonlinear equation (j52p . satisfied by the 
reduced one-body distribution function. (The latter in- 
dicates that a necessary condition to establish the ir- 
reversibility is the introduction of the thermodynamic 
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limit N —^ +00 .) The solution to Eq. (|52|) . in turn, de- 
termines the hierarchy of correlations namely Eq. ([70]) . 
We remark that although Eq. ((52|) is exact ^Jt differs 
from the usual relativistic kinetic equations [lol . [27j in 
the non-Markovian feature rendering the failure of H- 
theorem j40| . This feature was known long time ago 
in the Newtonian physics (fo r example, the Prigogine- 
Resibois master equation |5l|). It may, as pointed out 
in a pioneering study [i^, play essential roles in under- 
standing the infrared divergence of the collision integral 
of self-gravitating systems (which we will discuss below) . 
It is only beyond the correlation scale (provided that 
the two-body correlation is short-ranged) that the gen- 
eral equation (f52|) converges to a usual kinetic equation. 
(The very latter leads to a macroscopic hydrodynamic 
description.) In Paper II we exemplify this scenario in 
the case of relativistic plasmas with electromagnetic in- 
teractions. There, it will become clearer that Eqs. ([5^ 
and (|70p allow one to overcome the difficulty encountered 
in the Israel-Kandrup formalism. Furthermore, we man- 
age to justify the manifestly covariant Boltzmann equa- 
tion [13, [m , and are able to explore-at the quantitative 
level-the origin of the logarithmic divergence of the rela- 
tivistic Landau equation. They both are extremely hard 
to be achieved within the Israel-Kandrup formalism. 

It should be stressed that the general theory, under 
Assumption 2.1. and 3.1, is rigorous. These two assump- 
tions essentially require that (classical) chaotic dynamics 
sufficiently develops for the underlying relativistic many- 
body system. Similar to their Newtonian analogues to 
justify them mathematically is an extremely hard task 
and is unnecessary at this stage. Rather, from the prac- 
tical viewpoint the theory presented here may be appli- 
cable to a large class of systems where the interacting 
force preserves the mass-shell constraint and, moreover, 
does not depend on the acceleration of the acted particle. 
A heuristic example is a system composed of identical 
particles which interact with each other through a weak 
enough massive scalar field. In this case the interacting 
force is 



A. Relativistic self-gravitating systems 

In this series of papers we limit ourselves to the spe- 
cial relativity. The theoretic scope presented here may 
be further extended to the general relativity and, thus, 
find important applications in astrophysics. In particu- 
lar, there are no difficulties which impede generalizing the 
present manifestly covariant classical correlation dynam- 
ics to self-gravitating systems composed of relativistic 
star clusters or galaxies interacting merely gravitation- 
ally. Since this subject is far beyond the scope of this 
series of papers we here report briefly some preliminary 
observations and leave the thorough analysis for future 
studies. 

In formulating the classical correlation dynamics of rel- 
ativistic self- gravitating systems an additional difficulty 
arises, namely to treat the "gravitational force" . We here 
follow the prescription of Israel and Kandrup [ll] . First 
of all, we fix a background geometry (Al^,g^^(x)) . The 
8-dimensional /i phase space [compared with Eq. ([7])] 
associated with particle, say i is now defined as fii : 
{{xi,pt)\xi e M'^ ,g^y{xi)p^p''i = w?} . The background 
metric Qiivix) is such chosen that it solves some field 
equation with smoothed-out matter distribution as the 
source. Crucially, given particle i the realistic path Xi[<i\- 
dictated by underlying classical many-body dynamics-is 
the geodesic generated by another metric g'^jj (x) . It gen- 
erally deviates from the geodesic in the background ge- 
ometry (tW^, g^y(a;)) (and, as such, do not preserve the 
phase volume of /i^ .) This very deviation is driven by 
the difference of the Christoffel symbols F^^^ and F'^j^ 
(associated with the metrics g^^ and g'^^ , respectively), 

i.e., 5T^j,{x) = F'^^(a;) — F^j,(a;) , and describes the grav- 
itational force completely. Provide that 5F^j^(x) varies 
over a scale much smaller than the radius of the space- 
time curvature, it may be determined by some linear field 
equation. 

To substantiate these observations we need to make the 
following replacement in Eq. (|2ip (A^^^ — g^i, — Uf^iUui): 



Ft,{x,,P,) = Ao iv'^" - «)a..$.,(a;)|^^^^ , (71) 
(a^a^ + K^) $,,(.t) = -47rAo / dTjS'^^\x - .T,(r,)) , 



N r 



where |Ao$ij| <C m, the interaction strength is oc Ag , 
and K is the mass of the scalar field. Notice that the 
interacting force preserves the mass-shell constraint, i.e., 
Ui ■ Fij{xi,pi) ~ 0, but does not satisfy the conserva- 
tive condition namely Eq. ([^0]) . Nevertheless it turns out 
that the general theory well applies in this case and, fur- 
thermore, the entire program of Paper II may be readily 
carried over to this system. In particular, because of the 
short-ranged nature of the interacting force a well defined 
Landau collision integral results under the weak coupling 
approximation. 



+ mT^^^{xi)uxiUi 



d 



dpu 



(72) 



\L[^ - m5TU^,J)J—^'^^u'(ut + (z ^ j) , (73) 

where in Eq. ([75)1 the second term is obtained by ex- 
changing the particle labels of the first term. Here we 
notice that because of the linear field approximation men- 
tioned above the perturbed Christoffel symbol of particle 
i namely 5T^j^{xi) consists of iV — 1 contributions from 
all the other particles, and (5F^^(i, j) is the one associated 
with particle j . 

Then, one may further proceed to formally carry out 
the entire program presented here. In particular, we ob- 
tain, similar to Eq. (|52p . an exact closed nonlinear equa- 
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tion satisfied by the one-body distribution function: 

d 



fi2 



d^P2 
-i 



AA2/(2)k(l)=Kg[/],(74) 



where / is normahzed according to 
\imN ^+00 J^dY^fi J d'^pu'^'' f{x,p) — 1 with E a 
spacelike 3-surface, and IKg[/] is the collision integral 
(The explicit form may be readily found which is yet 
unnecessary for present discussions.). Because the 
perturbed Christoffel symbol is weak enough to the 
leading order A-expansion Eq. (j74p then gives 



<(9pi +mr^^(xi 



d 



dp^ 



'Ml 



/(1) = 0, (75) 



the solution to which self-consistently determines the 
background geometry, i.e., (A^^, giiu{x)) through the Ein- 
stein field equation [G^u is a function of the Christoffel 
symbol.): 



G^ = -8^ j ^u^p,fix,p) 



(76) 



Eqs. ([75]) and ([75)1 constitute the well-known Boltzmann- 
Einstein equations [1^]. 

Let us further replace X(^x,p) [cf. Eq. (|Tf)) ] by the phase 
trajectory generated by the geodesic in the background 
geometry. Provided that the self-gravitating system is 
dilute, by taking into account the so-called locality as- 
sumption we are able to show that the collision integral 
is 



Af] = 





(— 


-—) 




[dpt 


dpU 



d _ 

'Xf{xi,Pl)f{xi,P2) , 



d 
W2 



(77) 



2m'^[l - 2{ui ■ U2)^]^ 

/d*k 
S(k ■ ui)6(k ■ U2) 
V9 



Jkltf 



up to the accuracy (namely the weak coupling approx- 
imation to be detailed in Pap er II). This is the Israel- 
Kandrup collision integral ^M, well-known to astrophysi- 
cists. 

The collision integral (ffT]) justifies that the classical 
correlation dynamics can be well applied to relativistic 
self-gravitating systerns, which was questioned by Kan- 
drup long time ago [1^. It should be stressed, how- 
ever, that such a justification is at the level of formal 
manipulations. The serious difficulty is that this colli- 
sion integral suffers from the infrared divergence, and 
may be amounted to the failure of the locality assump- 
tion. (The ultraviolet divergence is due to the failure 
of the weak coupling approximation which may be read- 
ily healed and, thus, of no special interests [ia|-) Such 
a divergence, like plasmas [3q|, finds its origin at the 



long-ranged nature of gravitational forces, which was no- 
ticed by Chandrasekhar long time ago |38|. However, 
in (Coulomb) plasmas the Debye screening renders the 
effective potential short-ranged. As a result, the two- 
body correlation is short-ranged justifying the locality 
assumption, which results in a well defined kinetic equa- 
tion (One of the central issue of Paper II, indeed, 
is the detailed analysis of the analogue in the manifestly 
covariant classical correlation dynamics.) This scenario 
nevertheless breaks down in self-gravitating systems. In- 
deed, experiences in Newtonian self-gravitating systems 
suggest that the screening is dynamical [4l| . The general 
formalism presented here may serve as a useful technique 
for exploring this fundamental issue, which we leave for 
future studies. 
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APPENDIX A: BBGKY HIERARCHY OF 
PHYSICAL DISTRIBUTION FUNCTIONS 

Let us introduce the concept of the physical distribu- 
tion vector as follows: 

^ = ({M},{AA2},--- ,{AAjv}=AA). (Al) 

The general component above is defined as 

A^s is^i-i ^ Pii ^ ' ' ' 1 "^is 1 P^s ' '^il [^] ' ' ' ' ' [^] ) 

N 

(A2) 



H I di:^,^d^p,X. A/" 



for \ <ii<---<is<N,l<s<N. The normaliza- 
tion of the component is given by 



xAfs{xi^,Pi^, - ■ ■ ,Xi,,pjja;iJ<;], • • • ,2;^^^]) = 1 (A3) 



for arbitrary spacelike 3-surface . 

Similar to Sec. IIII Al from the Liouville equation 
we obtain 



(A4) 



Upon projecting it to the general s-particle state 
(ii , • • • , ?s| we use the same matrix elements except that 
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FIG. 1: (a) and (b): two types of interaction vertex, where 
dotted line stands for the annihilated particle, (c)-(f): exam- 
ples of diagrams with s' = 3 . 

the annihilation vertex (see Appendix [B|) is replaced by 
(«i , • • • ,«j|A£'|ii , • • • (A5) 

fe=l 

Notice that the BBGKY hierarchy (|A4p is manifestly co- 
variant. If or does not depend on the particle world lines, 
then we recover a nonmanifestly covariant BBGKY hier- 
archy by enforcing ti = ■ ■ ■ = tN = t , which is perfectly 
legitimate by relativity principles. 



\C[j ,i ^ j either preserves the particles (Fig. [T] (a)) or 
annihilates one particle via the manipulation of integra- 
tion (Fig.[T](b)). That is, 

(zj|A4-|r2(z,j)> - a4, 
(i|A4.|r2(z,j)) - j d]\c[^, 

(j|A4-|r2(z,j)) = j di\t[^, (B2) 

and all the other matrix elements vanish. In gen- 
eral, using the linearity we find that for matrix el- 
ements (rs'(ii,--- , Zs')|A£'|rs(ji, • • • not to van- 
ish it is necessary that (ji,--- ,js) = {ii,''' j*s') or 
(ii, • • • , is' ,is'+i) ■ Besides, provided that the correlation 
pattern |r,(ji, • • • , j,)> = |Pi|P2| • • • |Pj> with 

PiU---UPj = (ji,j2,--- ,is), 

P,nP,. = 0,Vi^i', (B3) 

then the matrix element is invariant under the particle 
permutation if both particles are in the same subset. This 
statement is also applicable to the correlation pattern 

Table 1 presents the matrix elements 
with s' = 1,2. For the matrix elements 

(r,,(ii,--- ,i,0|A£'|r,(ji,--- with s' > 2, we 

set up the diagrammatical rules as follows: 

1. The correlation pattern |rs(ji, • • • , js)) is given by 
the right-most part of the diagram. If particles are 
correlated then we connect them. 

2. For each particle draw a horizontal line which 
stands for the propagation from the right to the 
left. If particle is annihilated during the propaga- 
tion it is drawn by a dotted line, otherwise by a 
solid line. 



APPENDIX B: LIOUVILLIAN IN 
CORRELATION PATTERN REPRESENTATION 

In this appendix we will give the matrix elements of the 
free and interacting Liouvillian in the correlation pattern 
representation. For this purpose it is sufficient to know 
the matrix elements of the operators —u^d^i and A£'y 
because of the linearity. Let us project the BBGKY hi- 
erarchy (1^ to the correlation pattern (rs(ii,--- ,is)\ 
then substitute the cluster expansion psp into it. Com- 
paring with Eq. (pij) (setting A = 0) we find that the 
operator —u^d^i preserves a given correlation pattern 
|rs(ii, • • • , is)) (namely the free motion). That is, 

{Vs{ii, ■ ■ ■ ,is)\ (-u'^d^i) \Tsiii, ■ ■ ■ ,is)) = -"f^^i 

(Bl) 

if i G (ii, • • • , is) . All the other matrix elements vanish. 

Similarly, we find that, given a two-particle correlation 
pattern |r2(i,j)), the two-body interacting Liouvillian 



3. Draw all diagrams with an interacting vertex which 
may be either Fig.[T] (a) or (b). The obtained dia- 
gram, as a whole, gives a possible final correlation 
pattern which is read out according to the follow- 
ing rule: If particles are connected (irrespective of 
solid/dotted line) then they are correlated. 

4. Compare each possible final correlation pattern 
with (Fg' («!,••• ,is')\- If they are not identical 
then we assign the value zero to the correspond- 
ing diagram. 

5. Otherwise, depending on the type of the interac- 
tion vertex we assign the value XC'^j or / dk XC'^j^ 
accordingly. Here i,j,k are the particles joining the 
vertex. 

6. Summing up all the nonvanishing diagrams then 
gives the matrix element. 
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Table 1. Matrix element (rs/|A£'|rs) with s' = 1,2. All the other matrix elements not listed here are zero. 





m 




{i\X£'\i) =0 
(i\X£.'\i\j) = fdj XL' ■ 
(z|A£'|zj) = /djA4. 


mxsi'\t\j) = 

{i\j\X2'\ij) = 
mX£'mk)^JdkXiCi,+C'^„) 
{z\j\X£'\m=fdkXC'^, 

{i\j\xii'\m^jdkx£'^, 
mx£'m) = o 

{i\j\X£'\ijk) = 


{tj\X2'\t\j) = A£^ 
(ij\X£,'\ij) = A£'- 

{ij\x2'mk) = 
{ij\x£'m)^ jdkxci, 

{ij\X£'\j\ik)=JdkXCr, 
(zj|A£'|fc|*j)=/&A(£^,+4,) 
{tj\X£'\tjk) = JdkX{C'^,+£r,) 



For illustrations here we give some examples: and (ijfc|Aii'|i|j|A:) = . 
mX£'\t\jk) = A(4. + C[,) (Fig. m (c) and (d)), 

{ijk\Xt'\ij\kl) = (Fig. [11(e) and (f)), 
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